LINGS  S02  STUDY 

nalysis.  Model  Valida  tion 

and 

urce  Determlnatioii 

(DRAFT) 


VOLUME  I 

Text 


Ji&TE  CQUiCBW 

•;  fKiya  51284 

lAONI'ANA  SV&fS* 

1515  6.  ait  AVE, 
IIEiSHAj'MOHT  At*  A ' 5&S20 


■ ■ ■*,-  'j " 


Air  Quality  Bureau 
Montana  Department  of  Health 
and  Environmental  Sciences 


April,  1984 


JAN  2 5 1382 


BILLINGS  S02  STUDY 

Data  Analysis,  Model  Validation 
and 

Source  Determination 
(DRAFT) 


VOLUME  I 
Text 


BY 

Stan  Sternberg 
Ai r Qual i ty  Bureau 

Montana  Department  of  Health  and  Environmental  Sciences 

Cogswell  Building 
Helena,  Montana  59620 


Apri 1 , 1984 


Digitized  by  the  Internet  Archive 

in  2015 


https://archive.org/details/billingssostudyd1984ster 


PREFACE 


This  report  is  divided  into  three  volumes  for  printing  and  distribution 
purposes.  Volume  I contains  the  text  of  the  report.  It  includes  narratives  on 
data  analysis,  model  validation  and  source  determinations. 

Volume  II  contains  tables  and  figures  that  are  referenced  in  the  text.  The 
lists  of  tables  and  figures  are  repeated  in  this  volume  for  easy  reference. 

Volume  III  is  an  appendix  to  the  text  which  contains  a display  of  the  air 
quality  and  meteorological  data  collected  during  the  study  year.  Supportive 
material  is  also  provided  in  this  separate  volume. 
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ABBREVIATIONS 


0800  061  — (B061 ) 
0080  062  — (B062) 
0080  063  — (B063) 
0080  064  — (B064) 
0080  065  --  (B065) 
0080  066  --  (B066) 
0840  014  --  ( L014) 
0840  015  — (L015) 
0070  062  -- 


--  Metra  (Fairgrounds) 
--  Taft  School 
--  Shawnee  Park 
--  North  Johnson 
---  Lockwood  Park 
--  Coburn  Road 
--  New  Farm 
--  laurel-BN 
--  Taft  Composite 


42401  ---  Sulfur  Dioxide  - ppm 

61101  --  Wind  Speed  - m/sec 

61102  --  Wind  Direction  - Degrees 

61103  --  Compass  Wind  Direction  - Alpha 
61106  --  Standard  Deviation  of  Horizontal 
62101  --  Temperatre  - °C 

61501  --  Stability  --  Pasquil 1 -Gifford 
61401  --  Mixing  Height  - Meters 


Wind  Direction 


ppm  --  Parts  Per  Million 
ppb  --  Parts  Per  Billion 
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ug/m  Micrograms  per  Cubic  Meter 


AQDHS  --  Air  Quality  Data  Handling  System 
EPA  --  Environmental  Protection  Agency 
JCL  --  Job  Control  Language 
Met  --  Meteorology/Meteorological 
NWS  --  National  Weather  Service 

SAR0AD  Storage  and  Retrieval  of  Aerometric  Data 

S0^  --  Sulfur  Dioxide 

SPSS  --  Statistical  Package  for  the  Social  Sciences 

UNAMAP  --  User's  Network  for  Applied  Modeling  of  Air  Pollution 

USGS  --  United  States  Geological  Survey 

UTM  --  Universal  Transverse  Mercator 
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I.  INTRODUCTION 


This  report  has  been  prepared  by  the  Montana  Department  of  Health  and 
Environmental  Sciences  for  the  six  industries  in  Billings  and  Laurel  (signees  of 
the  1977  stipulation)  and  for  the  public  at  large. 

The  report  presents  an  analysis  of  the  Billings  SO^  and  meteorological  data 
collected  from  October,  1981  to  September,  1982.  A variety  of  means  are  used  to 
undserstand  the  data  and  present  it  in  summary  and  graphic  fashions. 

During  the  study  period,  the  six  industries  supplied  emission  data  to  the 
Department.  With  this  emission  data  combined  with  the  air  quality  and 
meteorological  data,  the  Department  attempted  to  validate  a dispersion  model. 
This  attempt  proved  unsuccessful.  Consequently,  the  Department  used  direction 
analysis  and  mass  emissions  to  show  source  contributions  by  the  Billings  and 
Laurel  industries  at  receptors  (monitors)  in  the  area. 


II.  BACKGROUND 


As  a result  of  a 1977  agreement  between  the  Montana  Department  of  Health 
and  Environmental  Sciences  and  industries  in  Billings  and  Laurel,  an  air 
monitoring  program  was  started  in  September,  1981.  Monitoring  was  conducted  at 
four  S02  sites  and  one  upper  air  site  by  Rockwell  International  from  September, 
1981  to  September,  1982.  The  Department's  Air  Quality  Bureau  also  monitored  S0;2 
at  four  other  sites  at  the  same  time.  Wind  speed  and  direction  data  was  also 
collected  at  the  SO^  sites.  Table  1 lists  the  air  monitoring  sites  and  Figures 
1 and  2 show  their  locations. 

The  industries  that  signed  the  stipulation  (or  agreement)  include  Montana 
Power  (coal-fired  power  plant);  Exxon,  Conoco  and  Cenex  (oil  refineries);  Great 
Western  Sugar  (sugar  beet  factory)  and  Montana  Sulphur  and  Chemical  (sulfur 
plant).  Cenex  is  located  in  Laurel  while  the  rest  are  in  Billings.  Figures  1 
and  2 show  the  locations  of  the  industries. 

The  stipulation  stated  that  the  monitoring  was  needed  to  determine  air 
quality  levels  of  SO^  (to  compare  against  standards).  The  stipulation  also 
stated  that  the  data  collected  (air  quality  and  meteorological)  would  be  used  by 
the  Department  to  validate  a model  and  determine  source  contributions  to  the 
ambient  levels. 

During  this  period,  the  industries  provided  monthly  reports  of  their 
emissions  to  the  Department.  Toward  the  end  of  the  study  some  stack  sampling 
was  conducted  by  a contractor  for  EPA. 
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III.  STUDY  YEAR 


Although  data  was  collected  from  September,  1981  to  September,  1982,  this 
report  only  dealt  with  the  time  period  from  October  1,  1981  to  September  30, 
1982.  Annual  averages  that  have  been  calculated  for  this  report  reflect  this 
time  period. 

The  Department  did  use  monthly  emission  data  that  was  reported  prior  to 
October,  1981  to  calculate  average  emissions. 
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IV.  FOUR  STUDY  MONTHS 


Due  to  the  mass  of  data  and  the  need  to  develop  a pattern  in  the  data,  the 
Department  decided  to  study  four  months  in  detail.  They  are: 

a.  January,  1982 

b.  May,  1982 

c.  August,  1982 

d.  November,  1981 

The  months  were  picked  for  seasonal  patterns  and  data  completeness. 
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V.  AMBIENT  SULFUR  DIOXIDE  DATA  SUMMARY 


The  first  approach  to  understanding  the  Billings  SO^  situation  was  to 
review  and  summarize  the  ambient  data.  The  data  is  stored  on  the  Department's 
Air  Quality  Data  Handling  System  (AQDHS)  and  can  be  retrieved  and  printed  out  in 
detail  lists.  These  lists  have  been  provided  to  the  industries  and  are  not 
included  in  this  report. 

Several  ambient  data  summaries  were  made  from  the  detail  lists  for  each  of 
the  eight  monitoring  sites.  Tables  2 through  13  provide  summaries  for  daily 
maximum  one-hour  readings.  Tables  14  through  25  provide  summaries  of  daily 
24-hour  average  concentrations.  Table  26  shows  monthly  SCL,  averages  and  Table 
27  shows  annual  averages. 
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VI.  DIURNAL  S02  PLOTS 

In  order  to  visualize  trends  in  the  SO^  data,  the  Department  made  diurnal 
S0?  plots.  Figures  3 through  22  show  diurnal  plots  for  most  of  the  sites  during 
the  four  study  months.  These  plots  show  each  individual  reading,  list  the 
highest  and  second  highest  reading  for  each  hour,  and  show  the  number  of 
observations  and  the  hourly  average.  The  hourly  averages  are  connected  with  a 
solid  line  to  show  the  diurnal  trend. 

In  general,  many  of  the  plots  show  higher  concentrations  during  midday 
(6-1400  Hours).  The  readings  also  appear  to  be  more  scattered  during  this 
period.  The  Coburn  Road  site  appears  to  stay  elevated  in  early  morning  but  then 
drops  down  in  the  afternoon.  The  patterns  seem  to  be  similar  (albeit  higher  or 
lower)  during  different  months. 
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VII.  EMISSIONS  FROM  INDUSTRIES 


Stack  data  and  SO^  emissions  were  supplied  to  the  Department  in  the  form  of 
monthly  reports.  An  agreement  was  made  before  the  study  for  the  industries  to 
only  report  data  for  the  second  Tuesday  of  each  month. 

Much  of  the  data  was  in  forms  that  could  not  be  used  as  direct  input  into  a 
model.  Most  problems  dealt  with  units  and  volume  corrections.  The  Department 
made  calculations  to  convert  this  data  into  a useable  form. 

Since  the  data  varied  from  month  to  month , averages  were  made  of  SO^ 
emissions,  flow  rates,  and  stack  temperatures.  It  was  felt  that  variations  from 
month  to  month  would  be  similar  to  variations  from  day  to  day  and  that  the 
averages  would  probably  be  representative  of  the  entire  year. 

Table  28  lists  the  average  emissions,  stack  data,  UTM  coordinates,  and 
elevations  in  metric  units.  Table  29  expresses  the  same  information  in  English 
units.  Figure  23  shows  a profile  of  stacks  and  receptors  with  elevations. 
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VIII.  STACK  SAMPLING  BY  GCA 


At  the  end  of  the  SO^  monitoring  period,  the  EPA  contracted  with  GCA  to  do 
stack  sampling  at  Exxon,  Conoco  and  Montana  Sulphur.  The  tests  were  performed 
during  October  and  November,  1982.  The  final  reports  were  received  by  the 
Department  in  December,  1983. 

It  was  hoped  that  this  emission  data  could  be  tied  with  that  from  the 
monthly  reports.  However,  after  review  by  the  Department,  problems  were  found 
with  the  flow  rates.  Consequently,  the  stack  tests  could  not  be  tied  in  with 
the  monthly  reports. 
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IX.  OPERATION  SCHEDULE 


The  Department  compiled  a listing  of  the  operations  of  21  sources 
(individual  stacks)  at  the  six  industries  in  Billings  and  Laurel.  Daily 
operation  schedules  are  presented  in  Tables  30  through  41.  Table  42  gives  an 
annual  summary. 

The  operation  schedules  were  used  in  writing  the  scenarios,  running  the 
models  and  calculating  mass  emissions. 

Montana  Power  and  Great  Western  Sugar  had  on/off  type  operations,  while  the 
oil  industries  and  Montana  Sulphur  ran  fairly  constantly  all  year. 
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X.  SELECTION  OF  MODELS 


Two  models  were  selected  to  try  to  validate  with  the  Billings  SO^ 
monitoring  data.  These  were  Complex  I and  Complex  II. 

Complex  I and  II  are  multiple  point  source  models  with  terrain  adjustments. 
They  are  sequential  models  utilizing  hourly  meteorological  input  (CRSMET  or 
RAMMET  type  preprocessor  data).  Complex  I assumes  that  hourly  averaged  plumes 
have  a normal  distribution  in  the  vertical  and  a uniform  distribution  across  a 
22.5  degree  sector  (similar  to  the  Valley  model).  Complex  II  assumes  the  plumes 
to  have  normal  distribution  in  both  the  horizontal  and  vertical. 

Complex  I and  II  are  the  only  models  in  the  EPA  UNAMAP  series  that  handle 
complex  terrain  and  take  hourly  sequential  data.  The  Billings  area  is  in 
complex  terrain  and  receptors  are  located  in  elevated  areas  as  well  as  the 
valley  floor.  It  is  possible  to  make  a preprocessor  tape  (hourly  sequential 
meteorological  data)  from  the  meteorological  data  collected  at  the  Billings  and 
Laruel  air  monitoring  sites. 
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XI.  PREPROCESSOR  MET  TAPE 


A preprocessor  met  tape  was  made  for  input  to  the  Complex  I and  II  models. 
The  customary  method  used  by  modelers  is  to  obtain  NWS  surface  observations  and 
twice  daily  mixing  heights  for  a whole  year  from  the  nearest  weather  station(s) 
and  to  process  that  data  using  the  CRSMET  program.  Output  from  this  program 
includes  year*  month,  day,  stability,  wind  speed,  flow  vector,  temperature  and 
mixing  height. 

The  Department  has  a program  similar  to  CRSMET  called  PRECRSTER.  PRECRSTER 
uses  site  specific  data  in  SAROAD  format  and  makes  a preprocessor  met  tape  like 
CRSMET.  PRECRSTER  requires  five  parameters  for  input.  These  are:  Wind  speed, 

wind  direction,  temperature,  stability  and  mixing  height.  Output  from  this 
program  is  the  same  as  from  CRSMET. 

The  modeling  programs  only  allow  one  wind  data  set  to  be  used.  Since  wind 
data  was  collected  at  eight  sites,  one  site  had  to  be  chosen  (or  sites 
composited)  to  be  representati ve  of  the  valley.  Before  picking  a site  the 
Department  decided  to  look  at  the  data  and  see  if  any  correlation  existed 
between  sites.  Table  43  shows  the  average  monthly  wind  speeds  at  the  eight  wind 
sites.  Wind  speeds  at  the  Coburn  Road  and  Laurel  sites  were  faster  than  at  the 
Billings/Lockwood  valley  floor  sites. 

The  Department  performed  regressions  of  wind  speeds  between  sites  and 
regressions  of  wind  directions  between  sites  for  two  random  months.  The 
correlation  coefficients  for  both  wind  speeds  and  wind  directions  were  between 
.5  and  .9.  Two  lobes  were  noted  on  the  scattergrams  for  wind  direction:  One  to 

the  NE  and  the  other  to  the  SW.  This  shows  up  and  down  valley  wind  flow 
patterns.  Wind  directions  between  sites  correlate  better  with  faster  winds 
(greater  than  2.5  m/sec)  than  with  slower  winds  (less  than  2.5  m/sec).  It  was 
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noted  that  wind  speeds  increase  with  increased  elevation:  Ridge  sites  have 

faster  winds  than  valley  floor  sites. 

Based  on  the  above  discussion,  the  Department  decided  that  wind  data  from 
one  site  would  be  representati ve  of  the  valley,  especially  if  long  term  modeling 
was  conducted  or  if  several  short  term  periods  were  averaged.  The  Department 
also  decided  to  pick  a valley  floor  site  rather  than  a ridge  site.  It  was  felt 
that  the  wind  profile  exponents  in  the  modeling  programs  would  take  care  of  the 
increase  in  wind  speeds  with  increased  elevation. 

The  Taft  School  wind  site  was  picked  as  the  representative  valley  floor 
site.  However,  as  is  the  nature  with  any  data  collection,  gaps  exist  in  the 
data  due  to  down  time  (repairs,  calibrations,  etc.).  Since  PRECRSTER.  requires  a 
complete  set  of  input  data,  it  was  necessary  to  fill  in  the  data  gaps  based  on  a 
heirarchy  of  sites.  The  hierarchy  is  as  follows: 


a. 

Taft  School 

- 

062 

b. 

Lockwood 

- 

065 

c. 

Metra 

- 

061 

d. 

Coburn 

- 

066 

e. 

Laurel  BN 

.. 

015 

Based  on  the  above  order,  gaps  were  filled  for  wind  speed,  wind  direction, 
and  temperature.  All  wind  speeds  less  than  one  m/sec  were  changed  to  one  m/sec. 
(In  modeling  formulas,  concentrations  are  proportional  to  the  inverse  of  wind 
speed.  If  wind  speeds  are  less  than  one  m/sec  or  zero,  the  formulas  calculate 
extremely  high  concentrations,  or  with  zero,  cause  division  problems.) 

Two  other  parameters  were  needed  for  input  to  the  preprocessor  program. 
These  were  stability  and  mixing  height.  PRECRSTER  requires  stability  to  be 
entered  in  two  hour  intervals  (12  per  day)  and  mixing  height  in  six  hour 
intervals  (four  per  day). 
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The  modeling  formulas  require  the  input  of  the  standard  Pasqui 11 -Gifford 
(PG)  type  stability  (A  through  F).  The  Department  reviewed  the  Rockwell  reports 
and  decided  that  the  stability  classifications  reported  could  not  be  converted 
to  a PG  classification.  Thus,  an  alternative  approach  was  employed.  The 
Department  used  the  CRSMET  program  and  Billings  airport  NWS  surface  observations 
for  the  study  year  to  determine  stabilities.  The  PG  stabilities  were  taken  off 
of  the  CRSMET  output  and  put  in  SAROAD  format  for  the  PRECRSTER  program.  The 
PRECRSTER  program  only  allows  single  increment  steps  of  stability  from  hour  to 
hour. 

The  Department  obtained  mixing  height  data  at  hours  of  0000,  0600,  1200  and 
1800  for  each  day  of  the  study  year  from  the  acoustic  sounder  data  presented  in 
the  Rockwell  report.  The  heights  were  then  put  in  SAROAD  format  and  used  in  the 
PRECRSTER  program.  PRECRSTER  reads  the  data  in  six  hour  intervals  and 
interpolates  for  hours  in  between. 

Rockwell  presented  two  default  values  or  height  codes  in  their  report.  The 
value  888  meant  unlimited  mixing  and  999  meant  missing  data.  The  Department 
filled  in  data  gaps  (999)  by  substituting  a mixing  height  that  occurred  close  to 
the  hour  that  was  missing,  or  using  a mixing  height  that  occurred  at  the  same 
time  on  a day  before  or  after  (assuming  similar  meteorology).  The  unlimited 
value  of  888  was  kept  at  888  meters.  The  Department  felt  that  this  mixing 
height  would  be  high  enough  so  as  not  to  effect  the  plumes  from  Billings 
industries  and  so  that  hourly  interpolations  would  be  fairly  realistic. 

In  review,  the  Department  made  data  sets  complete  for  one  year  for  wind 
speed,  wind  direction,  temperature,  stability  and  mixing  height.  The  program 
PRECRSTER  processed  the  data  sets  and  made  a tape  which  is  input  to  the  models. 
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An  example  of  input  to  PRECRSTER  is  shown  in  Figure  24  and  output  is  shown 

fc 

in  Figure  25.  An  example  of  input  to  CRSMET  is  shown  in  Figure  26  and  output  is 

shown  in  Figure  27. 
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XII.  STATISTICAL  ANALYSIS 


In  order  to  gain  understanding  of  the  Billings  situation,  the  Department 
applied  statistical  analyses  to  the  data  base.  The  Department  used  the  state's 
statistical  package  (SPSS)  for  this  analysis.  Before  application  of  SPSS,  the 
hourly  records  which  are  stored  in  horizontal  rows  (12  per  line  SAROAD  format) 
had  to  be  rearranged  into  vertical  columns.  The  Department  used  a program 
called  PRESTAT  to  do  this. 

Data  sets  were  made  for  each  of  the  four  study  months.  Each  data  set  had 
36  variables  by  744  or  720  hours  in  a month.  The  data  sets  contained  a counter, 
Julian  day,  month  and  day,  hour  of  day  and  then  S09,  wind  speed,  wind  direction, 
and  standard  deviation  of  horizontal  wind  direction  for  each  of  the  eight  sites. 
Figure  28  gives  an  example  of  part  of  one  data  set.  Using  these  data  sets 
statistical  analyses  (regression,  selection  criteria,  etc.)  can  be  run. 

Tables  44  to  47  show  means  and  correlation  coefficients  comparing  hourly 
SO^  data  between  sites  for  the  four  study  months.  There  appears  to  be  some 
correlation  between  sites  in  the  Lockwood  area  and  between  sites  in  Laurel. 

SPSS  and  PRESTAT  data  sets  were  used  to  make  data  sets  for  input  to  the 
direction  analysis  program.  Results  of  this  are  presented  in  the  next  section. 
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XIII.  DIRECTION  ANALYSIS 


The  direction  analysis  program  (DIRECTAN)  was  set  up  by  the  Department  to 
quantify  and  summarize  SO^  concentrations  monitored  at  specific  sites  based  on 
wind  direction  criteria.  Input  to  the  DIRECTAN  program  was  based  on  selection 
criteria  set  up  with  SPSS  and  using  the  PRESTAT  data  sets. 

An  example  of  the  selection  criteria  (in  SPSS  code)  is  as  follows: 

SELECT  IF  (BQ63WD  GT  101.25  AND  B063WD  LT  123.75) 

CONDESCRIPTIVE  B063S2  BQ63WS 

This  reads:  If  an  hourly  reading  of  wind  direction  from  the  Shawnee  Park 

site  (B063WD)  is  greater  than  101.25  degrees  and  less  than  123.75  degrees  (or 
from  the  ESE  compass  direction),  then  condescript  the  Shawnee  Park  SO^  readings 
(B063S2)  and  condescript  the  Shawnee  Park  wind  speed  readings  (B063WS). 

The  above  selection  criteria  was  done  in  turn  for  each  of  the  16  compass 
directions  (see  Figure  29)  and  for  each  of  the  sites  for  the  four  study  months. 
Condescript  reports  the  mean,  the  number  of  valid  observations  and  the  number  of 
missing  observations  for  SO^  and  wind  speed.  This  data  was  then  input  to  the 
DIRECTAN  program. 

The  DIRECTAN  program  displays  the  input  data  and  makes  calculations. 

Values  that  are  calculated  include  percent  wind  direction,  SO^  strength  factor 

and  percent  strength  factor. 

SO^  strength  factor  is  the  mean  of  the  S0,>  readings  (from  one  of  the  16 
compass  directions)  times  the  percent  wind  direction.  The  strength  factor  gives 
a weighted  SO^  concentration.  Percent  strength  factor  is  obtained  by  adding  up 
individual  strength  factors  and  then  dividing  each  by  the  total  times  100. 
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The  percent  strength  factors  show  the  contribution  from  each  compass 
direction  on  the  total  monthly  average.  Theoretically,  the  percents  can  be 
matched  to  the  sources  in  those  directions  and  source  contribution  shown. 

Tables  48  to  75  show  the  direction  analysis  results  for  sites  during  the 
four  study  months.  General  comments  are  made  in  regard  to  the  tables: 

a.  Sites  in  Laurel  show  two  lobes:  One  from  the  SW  and  the  other  from 

the  NE. 

b.  The  Taft  School  and  Coburn  Road  sites  show  two  lobes:  One  from  the  SW 

and  the  other  from  the  NE. 

c.  The  Lockwood  Park  site  shows  two  lobes:  One  from  the  north  and  the 

other  from  the  SW.  The  lobe  from  the  north  does  not  occur  too  often 
but  the  results  show  high  concentrations. 

d.  The  sites  to  the  east  of  Billings  (Shawnee  Park  and  North  Johnson) 
show  only  one  lobe  from  the  SW. 

e.  SO^  is  seen  coming  from  all  directions  at  all  sites.  This  might  be 
attributed  to  calm  conditions  where  direction  has  little  meaning  or 
due  to  inversion  breakup  times  when  SO^  is  trapped  above  at  night  and 
then  is  transported  down  when  the  sun  begins  to  heat  the  ground 
(independent  of  wind  direction). 
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XIV.  DISPLAY  OF  AIR  QUALITY  AND  METEOROLOGICAL  DATA 


In  order  to  review  the  data,  the  Department  made  a display  of  the  meteoro- 
logical data  with  the  ambient  S02  data  for  each  day  of  the  study  year.  The 
Department  thought  that  a display  of  this  sort  would  allow  a reviewer  to  look  at 
a day's  data  and  see  if  any  trends  existed  (visual  observation)  and  get  an  idea 

from  where  the  S02  was  coming. 

The  meteorological  data  that  was  used  for  this  display  came  from  the  met 
data  that  was  assembled  for  the  preprocessor  tape  (see  Section  XI).  The  wind 
speeds  that  were  less  than  one  m/sec  were  kept  as  originally  measured  and  not 
rounded  up  to  one  m/sec  as  was  done  for  the  preprocessor.  In  addition  to  the 
five  parameters  used  for  the  preprocessor  (wind  speed,  wind  direction,  tempera- 
ture, stability  and  mixing  height),  standard  deviation  of  horizontal  wind 
direction  and  the  alpha  characters  for  the  16  compass  directions  were  added  to 
the  display.  As  a naming  convention,  the  Department  refers  to  this  met  data  set 
as  Taft  Composite. 

The  S02  data  for  each  of  the  eight  sites  was  retrieved  from  the  Depart- 
ment's computer  system  and  sorted  by  date,  site  and  parameter  with  the  met  data. 
The  display  uses  the  Department's  minimaster  report  format  which  presents 
labels,  spreads  out  the  SAROAD  record,  and  inserts  a decimal  for  ease  of  read- 
ing. 

Highlights  were  added  by  hand.  Brackets  were  added  around  each  site's  24 
hour  readings.  A number  was  added  to  the  left  of  the  brackets  which  gives  the 
24-hour  average  in  parts  per  billion.  The  maximum  one-hour  readings  are 

underlined  for  each  site. 

The  displays  for  the  whole  study  year  are  presented  in  the  appendix.  The 
year  for  October,  November  and  December  is  1981  rather  than  1982  as  is  shown. 
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The  Department  reviewed  four  months  of  the  displays  and  presented  the  results  in 
the  form  of  daily  scenarios.  Results  are  presented  in  the  next  section. 
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XV.  SCENARIOS 


The  display  of  met  and  ambient  SO^  data  was  set  up  so  that  the  Department 
could  visually  look  at  the  data  and  write  a brief  scenario  of  what  was  occurring 
for  each  day  of  the  four  study  months.  A check  sheet  was  set  up  to  note  meteo- 
rology, Billings  area  SO^  trends.  Laurel  area  SO^  trends,  operation  schedules 

and  comments. 

For  S0o  concentration  trends,  the  Department  noted  hours  of  similar 
elevated  readings  or  peaks  that  stood  out  over  a background  level  or  baseline. 
For  the  same  time  period,  the  Department  noted  the  corresponding  meteorology  and 
the  subareas  or  sources  that  caused  or  contributed  to  the  concentrations.  There 
was  no  attempt  to  note  how  strong  the  elevated  readings  were,  but  just  their 
occurrence. 

The  letters  that  were  used  for  the  subareas  are  as  follows: 

A - Cenex 

B - Montana  Power,  Conoco  and  Great  Western 
C - Exxon,  Montana  Sulphur 

Tables  76  to  91  give  the  daily  scenarios  for  the  four  study  months.  Tables 
92  to  95  summarize  some  of  the  comments  from  the  daily  scenarios. 

The  following  points  are  made: 

1.  Inversion  breakup  and  buildups  during  calm  periods  are  major  factors 

that  cause  elevated  SO^  readings  to  occur  at  the  Billings  monitors. 

Breakups  cause  elevated  readings  to  occur  in  the  early  morning  to  mid-day 

(6-1400  hours). 

2.  SO^  can  be  seen  going  from  Billings  industries  to  Laurel  monitors  and 

from  Laurel  industries  to  Billings  monitors.  In  the  scenario  summaries, 

there  are  more  days  noted  when  SO^  traveled  from  Billings  to  Laurel.  This 
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does  not  mean  that  the  phenomenon  happens  more  than  the  reverse  but  that  it 
was  easier  to  note. 

3.  During  calm  periods  and  high  mixing  heights  (300-400  meters),  SC^  from 
the  stacks  rises  and  buildups  occur  aloft.  During  this  time  concentrations 
are  not  seen  at  the  valley  floor  monitoring  sites,  but  can  be  seen  at  the 
higher  Coburn  Road  site. 

4.  During  calm  periods  and  low  mixing  heights  (30-200  meters),  SO.,  from 
the  stacks  is  not  allowed  to  rise  too  far.  Buildups  can  occur  at  ground 
level  where  the  monitors  are  located. 

5.  Faster  wind  speeds  (greater  than  2-3  m/sec)  cause  the  plumes  from  the 
industrial  sources  to  impact  on  the  ground.  A combination  of  faster  winds 
dominating  plume  rise  and  blowing  the  plumes  parallel  to  the  ground  and 
terrain  effects  (causing  downwash)  may  be  the  reason  why  the  plumes  are 
coming  down. 

6.  When  there  are  faster  wind  speeds,  the  wind  direction  data  becomes 
more  reliable.  By  having  more  confidence  in  the  wind  direction  data,  it  is 
easier  to  tell  from  which  subarea  or  source  the  SO^  was  coming. 

7.  When  there  are  slower  winds,  the  wind  direction  data  became  less 
reliable  and  consequently  it  became  harder  to  tell  from  which  subarea  or 
source  the  S0^  was  coming. 

8.  Inversion  breakups  and  faster  winds  bring  the  plumes  down  and  can 
occur  at  the  same  time  and  have  an  additive  effect.  For  example,  buildup 
can  occur  at  night  and,  in  the  morning,  inversion  breakup  can  occur.  At 
the  same  time  in  the  morning  the  wind  speeds  start  to  pick  up  and  bring  the 
"just  emitted"  SO^  to  the  ground.  Together  the  old  emissions  that  were 
build  up  at  night  add  up  with  the  new  emissions  during  the  day  and  cause 
higher  concentrations  at  a site. 
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9.  Each  of  the  subareas  can  have  an  additive  effect  on  each  other  as  the 
wind  flows  either  up  or  down  the  valley.  For  example  with  a SW  wind,  S(L 
emitted  from  subarea  A can  first  be  seen  at  the  Laurel-BN  site.  Then  the 
SO^  travels  on  to  Billings  where  it  adds  up  with  emissions  from  subarea  B 
and  together  they  impact  the  Coburn  Road  and  Lockwood  sites.  Then  the  SC^ 
from  subareas  A and  B add  up  with  emissions  from  subarea  C and  are  seen  at 
the  North  Johnson  and  Shawnee  Park  sites. 

The  reverse  can  also  be  seen:  SC^  from  subarea  C can  be  emitted  with 

a NE  wind  and  is  measured  at  the  Coburn  Road  site.  Then  the  SO^  adds  up 
with  emissions  from  subarea  B and  both  go  on  to  Laurel  where  they  are 
monitored  at  the  Laurel-BN  and  New  Farm  sites.  It  would  be  assumed  that  SW 
of  Laurel,  subareas  C and  B would  add  up  with  emissions  from  A,  however,  no 
monitor  exists. 
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XVI.  DIURNAL  TREND  PLOTS 


In  order  to  gain  more  understanding  into  the  diurnal  trends  in  Billings, 
the  Department  made  diurnal  plots  of  SQ^  wind  speed  and  standard  deviation  of 
horizontal  wind  direction  for  each  site  for  the  four  study  months. 

Figures  30  to  57  show  these  plots. 

It  is  interesting  to  note  that  the  traces  for  S0^  and  wind  speed  follow 
each  other.  In  fact,  the  wind  speed  actually  lags  a little  behind  the  SO^ 
trace.  This  means  that  there  is  another  factor  that  causes  each  to  go  up.  This 
can  be  attributed  to  the  heating  of  the  ground  by  the  sun.  The  sun  comes  up  and 
heats  the  ground  which  causes  breakup  and  increased  wind  speeds.  In  general, 
sunrise  corresponds  to  the  time  of  breakup. 

When  the  sun  sets,  wind  speeds  die  down  and  concentrations  go  down. 

Buildups  then  occur  later  (20-2400  hours)  either  aloft  or  near  the  ground. 
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XVII.  WIND  ROSES 


The  Department  made  a series  of  wind  rose  plots  using  the  Taft  Composite 
wind  data  set.  These  wind  roses  are  an  aid  in  understanding  wind  flow  patterns 
in  Billings  and  along  the  Yellowstone  River  valley. 

Figure  58  gives  a wind  rose  for  the  whole  study  year.  Figure  59  gives 
quarterly  wind  roses.  Figures  60  to  62  give  monthly  wind  roses.  Figure  63 
gives  monthly  wind  roses  for  each  of  the  four  study  months.  Figures  64  to  71 
give  diurnal  wind  roses  for  the  four  study  months.  Tables  96  to  144  give  the 
wind  rose  frequencies  for  each  of  the  wind  roses. 

A few  comments  are  made  concerning  the  diurnal  wind  roses.  In  the  early 
morning,  the  winds  are  mostly  the  slow  river  drainage  southwest  types.  In  the 
afternoons  the  wind  directions  often  change  and  the  speeds  pick  up.  Then  in 
late  evening  the  winds  die  down  again  and  follow  the  river  drainage. 

The  monthly,  quarterly  and  annual  wind  roses  all  show  similar  patterns. 
There  are  essentially  two  nodes:  One  from  the  SW  and  the  other  from  the  NE. 

The  node  from  the  SW  occurs  about  twice  as  often  as  the  node  from  the  NE. 

Winds  from  the  SE  show  low  wind  speeds  and  reflect  calm,  variable  con- 
ditions. Winds  from  the  NW  usually  show  fast  wind  speeds  and  are  probably 
caused  by  frontal  passages. 
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XVIII.  MASS  EMISSION  CALCULATIONS 


Mass  emissions  (tons/month  or  tons/year)  are  calculated  by  multiplying  an 
emission  rate  (g/sec)  times  an  operation  period  (days/month  or  days/year)  and 
using  appropriate  conversion  factors.  Total  mass  emissions  in  an  area  (Billings 
or  Laurel)  can  be  found  by  adding  up  the  individual  mass  emissions.  Percents 
can  be  found  by  dividing  the  individual  emissions  by  the  total  times  100. 

The  emission  rates  were  presented  in  Section  VII  and  the  operation  periods 
were  presented  in  Section  IX. 

Tables  145  and  146  give  mass  emissions  by  ton  and  percent  for  each  source. 
Tables  147  and  148  give  mass  emissions  by  ton  and  percent  for  each  industry. 
Emissions  from  sources  in  Billings  and  Laurel  were  kept  separate  for  totaling 
purposes . 

Table  149  presents  the  results  of  a regression  that  was  done  between  the 
monthly  mass  emissions  from  the  Billings  sources  and  the  monthly  SO^  averages 
for  each  of  the  eight  monitoring  sites.  This  table  shows  correlation  coeffic- 
ients between  .6  and  .9  when  all  twelve  months  are  present.  Although  the 
results  show  good  correlation  between  emissions  and  ambient  concentrations  (over 
a long  term  period),  it  shouldn't  necessarily  be  assumed  that  mass  emission  is 
the  only  factor  effecting  concentrations  (especially  for  short  term  concentra- 
tions). 

The  effects  that  meteorology  and  terrain  cause  on  the  dispersion  of  plumes 
and  the  relative  positions  of  the  monitors  with  respect  to  the  industries  also 
play  important  roles  on  the  concentrations  that  are  seen  at  the  monitors  (both 
for  the  long  term  and  for  the  short  term). 

Mass  emission  calculations  can  (in  a rough  manner)  be  used  to  approximate 
source  contribution. 
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XIX.  MODELING  RESULTS 


As  mentioned,  the  intent  of  this  report  was  to  validate  a model (s)  and  to 
show  source  contribution  by  several  Billings  and  Laurel  industries  at  receptors 
in  the  area.  This  section  shows  what  was  done  to  set  up  and  run  the  models 
(options  and  input  data)  and  what  resulted  as  the  programmed  output.  Discus- 
sions will  center  around  what  was  done  to  validate  the  model (s)  and  why  they 
failed  to  match  up  with  the  ambient  data  collected  at  the  monitoring  sites. 

In  review,  the  Department  chose  the  Complex  I and  II  models  because  of 
Billings'  complex  terrain  and  set  up  a representative  preprocessor  met  tape 
which  would  be  input  to  the  models.  The  Department  used  stack  emissions  that 
were  supplied  to  it  by  the  industries  and  located  UTM  coordinates  and  elevations 
for  the  receptors  and  sources  on  USGS  maps  (see  Table  28).  The  operation 
schedule  presented  the  Department  with  the  knowledge  of  which  sources  ran  on 
which  days. 

The  four  study  months  were  selected  to  be  modeled.  The  Complex  I and  II 
models  were  run  on  identical  input  data  and  the  results  of  both  were  compared 

with  the  actual  ambient  data  as  well  as  each  other. 

It  was  desirable  to  run  the  models  for  one  month  at  a time  using  the  same 

set  of  stack  information  for  the  whole  month.  However,  since  some  sources  ran 

less  than  a full  month,  rounding  was  done  either  up  to  the  full  number  of  days 
in  the  month  or  down  to  zero  days.  The  following  was  done: 

a.  For  January,  the  Exxon  Cat  unit  and  the  Cenex  #3  Boiler  were  rounded 
up  to  31  days. 

b.  For  May,  the  Exxon  Cat  unit  was  rounded  to  zero  days.  Cenex' s sources 
were  rounded  to  31  days. 
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Co  For  August,  the  Montana  Power  Corette  plant  was  rounded  to  zero  days. 

Cenex 's  sources  were  rounded  to  31  days. 

d.  For  November,  no  changes  were  made.  All  sources  were  30  days. 

In  addition  to  the  point  source  information,  the  Department  selected 
several  options  and  parameters  for  input  to  the  models.  These  included: 

Starting  Julian  day,  starting  hour,  number  of  averaging  periods,  number  of  hours 
in  an  averaging  period,  technical  and  control  options,  and  wind  and  terrain 
factors. 

The  starting  Julian  days  for  each  month  were  as  follows: 


a. 

January 

1 

b. 

May 

121 

c. 

August 

213 

d. 

November 

305 

The  starting  hour  for  each  month  was  set  at  one;  the  number  of  averaging 
periods  corresponded  to  the  number  of  days  in  the  month  (30  or  31);  and  the 
number  of  hours  in  an  averaging  period  was  set  at  24.  (The  above  allowed  the 
models  to  calculate  twenty-four  hour  concentrations  for  each  site  for  each  day 
of  the  four  study  months.) 

The  Department  also  ran  the  Complex  II  model  to  calculate  a monthly  average 
for  January  and  to  show  modeled  source  contribution.  For  this,  the  starting 
Julian  day  was  one,  the  starting  hour  was  one,  the  number  of  averaging  periods 
was  one,  and  the  number  of  hours  in  an  averaging  period  was  744. 

The  technical  and  control  options  included  the  following: 

I OPT ( 1 ) = 1 Use  terrain  adjustments 

I0PT ( 2 )=1  Do  not  include  stack  downwash  calculations 

I OPT ( 3 ) = 1 Do  not  include  gradual  plume  rise  calculations 

I OPT ( 4 ) = 1 Include  buoyancy  induced  dispersion 
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I0PT(17)=0  Provide  averaging-period  contributions 

I0PT( 18)=0  Provide  averaging  period  summary 

IOPT ( 24 )= 1 Punch  averaging  period  concentrations  on  cards  (disk) 

IOPT ( 25 )= 1 Use  complex  terrain  option 

IOPT ( 24 ) = 1 allowed  averaging  period  concentrations  to  be  written  to  disk  in 
an  eighty  character  card  format.  The  actual  twenty-four  hour  concentrations 
were  then  laid  over  (or  next  to)  the  modeled  concentrations.  The  new  data  set 
was  then  used  in  SPSS  to  run  regressions  of  actual  versus  modeled  concentra- 
tions. 

Wind  and  terrain  factors  were  selected  for  input  to  the  models.  The 
Department  chose  wind  profile  power  law  exponents  that  reflected  the  higher  wind 
speeds  that  were  seen  at  the  Coburn  Road  site  versus  what  was  seen  at  the  valley 
floor  sites.  The  wind  profile  exponents  that  were  input  to  the  models  were  as 
follows:  .1,  .1,  .1,  .1.,  .1,  .1  for  each  of  the  stabilities  (A  through  F). 

The  terrain  adjustment  factors  selected  for  the  models  were  the  EPA  default 
values.  These  were  as  follows:  .5,  .5,  .5,  .5,  .0,  .0  for  each  of  the  stabil- 
ities (A  through  F) . 

Figures  72  to  75  show  the  input  data  sets  that  were  used  to  model  the  point 
sources  for  the  four  study  months.  The  figures  also  show  the  associated  JCL 
(Job  Control  Language)  that  goes  along  with  running  the  models. 

The  models  were  run  on  the  State's  IBM  3033N  computer  and  the  results 
consisted  of  printed  output  and  disk  data  sets  containing  the  averaging  period 
concentrations  for  each  site  and  for  each  day  of  the  four  study  months. 

The  printed  output  consists  of  summary  and  descriptive  information  (manda- 
tory output)  as  well  as  the  averaging  period  contributions  and  the  averaging 
period  summary.  The  mandatory  output  consists  of  the  following:  General  input, 
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point  source  data,  significant  point  source  ranking,  receptor  data  and  addi- 
tional information. 

Information  is  provided  in  the  mandatory  output  which  helps  to  focus  in  on 
how  the  models  handle  (and  apportion)  the  various  point  sources.  Table  150 
shows  the  significant  point  source  ranking,  chi-max  and  the  effective  height. 
Chi-max  is  the  maximum  expected  ground-level  impact  under  B stability  with  a 
stack  top  wind  speed  of  3 m/sec.  The  significant  point  source  ranking  is  based 
on  chi-max.  The  effective  height  is  the  sum  of  the  stack  height  plus  plume  rise 
(based  on  Briggs  plume  rise  formulas).  Generally,  sources  that  show  a large 
chi -max  (low  ranking  number)  are  usually  the  ones  that  the  models  say  are  the 
major  contributors.  Sources  with  large  effective  heights  are  usually  not  major 
sources  (as  treated  by  the  models)  because  the  plumes  are  allowed  to  go  up  and 
disperse  Into  the  atmosphere.  Later  discussions  will  focus  on  this 
information. 

The  averaging  period  summary  and  the  concentrations  written  to  disk 
(I0PT(24)=1,  one  data  set  each  for  Complex  I and  II)  gave  the  modeled  twenty- 
four  hour  concentrations  for  each  site  for  each  day.  Using  a computer  overlay 
function  and  the  I OPT ( 24 ) data  sets,  new  data  sets  were  made  (one  for  each  of 
the  study  months)  that  contained  the  Julian  day,  actual  concentration,  Complex  I 
modeled  concentration,  Complex  II  modeled  concentration,  and  site. 

Using  the  above  new  data  sets  and  SPSS,  regressions  were  made  between  the 
actual  and  modeled  concentrations.  Tables  151  to  154  show  the  means,  number  of 
pairs  and  the  correlation  coefficients  for  all  data  and  for  each  site.  General- 
ly the  models  underpredicted  (compared  to  actual  concentrations)  and  the  corre- 
lations were  poor. 

The  printed  output  also  gives  the  averaging  period  contributions.  The 
averaging  period  contributions  show  the  micrograms  per  cubic  meter  that  each 
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point  source  contributes  to  the  total  averaging  period.  When  the  averaging 
period  is  twenty-four  hours,  the  output  gives  averaging  period  contributions  for 
each  day  and  site.  However,  it  is  impractical  to  summarize  so  many  averaging 
period  contributions.  (There  are  21  sources  x 123  days  x 8 sites  = 20664  daily 

source  contributions.) 

Consequently,  for  demonstration  purposes,  the  Department  ran  the  Complex  II 
model  on  a monthly  basis  (mentioned  before).  Table  155  shows  monthly  modeled 
and  actual  averages  and  the  percent  contributions  from  the  sources  for  the  month 

of  January. 

The  percent  contributions  in  Table  155  are  not  realistic.  According  to  the 
model,  the  Montana  Sulphur  stack  seems  to  be  the  major  source  in  Billings.  In 
Laurel,  the  major  source  appears  to  be  the  Cenex  sulfur  plant.  Also,  there 
appears  to  be  hardly  any  contribution  by  the  Montana  Power  Corette  plant.  (This 
seems  unrealistic  since  the  Department  has  seen  plume  impact  by  this  source 
during  January  at  Coburn  Road  and  other  sites.)  The  percent  contributions  for 
the  monthly  averaging  period  are  similar  to  those  seen  for  the  shorter  twenty- 
four  hour  averaging  periods.  The  Department  concludes  that  the  modeled  percent 
contributions  can  not  be  used  for  source  apportionment. 

The  Department  tried  to  understand  why  the  models  (as  used)  were  unable  to 
match  up  with  the  actual  concentrations.  A number  of  parameters  were  looked  at 
which  effect  modeled  concentrations  and  also  how  they  effect  or  coincide  with 
concentrations  that  are  seen  in  the  field.  They  are  as  follows:  Effective 

plume  rise,  wind  speed,  fumigation,  and  receptor  location. 

The  effective  plume  rise  is  calculated  in  the  models  using  the  Briggs  plume 
rise  formulas.  For  example,  the  model  calculates  an  effective  plume  height  of 
485  meters  for  the  Montana  Power  Corette  plant  (stack  height  plus  plume  rise). 

In  reality,  the  plume  might  rise  from  30  to  60  meters  above  stack  height  (or 
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have  an  effective  plume  height  of  140  to  170  meters).  The  other  sources  in 
Billings  and  Laurel  may  show  patterns  similar  to  the  Montana  Power  Corette 
plant:  Actual  plume  rise  is  less  than  calculated  by  the  models  (Briggs 

formulas).  The  effect  that  the  plume  rise  has  on  modeled  concentrations  is  that 
sources  with  large  plume  rises  are  allowed  to  go  up  higher  and  disperse  (and  do 
not  contribute  to  modeled  concentrations  at  the  receptors). 

In  modeling  formulas,  concentrations  are  inversely  proportional  to  wind 
speed.  The  logic  being:  The  faster  the  wind  blows,  the  more  mixing  and 

dispersion  takes  place  and  lower  concentrations  are  seen  at  the  receptors. 
However,  in  Billings  quite  different  situations  take  place.  With  calm  winds, 
the  emissions  from  the  stacks  rise  and  concentrations  are  not  seen  at  the 
receptors  (unless  concentrations  build  up  over  the  whole  valley  with  a low 
mixing  height  and  then  they  are  seen  at  the  receptors). 

With  faster  wind  speeds  (greater  than  2-3  m/sec),  the  plumes  from  the 
sources  come  down  and  impact  the  receptors  (direction  dependent).  This  was 
observed  by  looking  at  the  met  and  ambient  data  for  the  scenarios. 

Another  case  exists  with  faster  wind  speeds.  When  inversion  breakup  occurs 
(due  to  ground  heating)  and  concentrations  increase,  wind  speeds  increase 
simultaneously  (due  to  the  same  ground  heating).  The  net  effect  being:  No 

relationship  between  SO^  concentrations  and  wind  speeds  occur  as  far  as  modeling 
is  concerned. 

Another  reason  why  the  models  are  unable  to  predict  accurate  concentrations 
in  Billings  is  because  of  their  inability  to  handle  fumigations  or  inversion 
breakups.  The  following  quote  from  the  MPTER  User  Guide  discusses  this  problem: 
"MPTER  (Complex  I and  II)  does  not  include  calculations  for  the 
transitional  phenomenon  of  fumigation,  which  is  the  elimination  of  an 
inversion  layer  containing  a stabilized  plume  from  below  (which  causes 
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mixing  of  pollutants  downward,  resulting  in  uniform  concentrations 
with  height  beneath  the  original  plume  centerline)." 

As  a result,  the  models  do  not  have  any  memory  of  what  has  been  emitted  and 
built  up  at  night  and  then  allowed  to  come  down  and  disperse  at  midday. 

The  locations  of  the  receptors  with  respect  to  the  sources  may  also  be  a 
reason  why  the  models  do  not  accurately  predict  concentrations.  The  receptors 
are  fairly  close  to  the  sources  and  the  model's  source-receptor  relationship  is 
highly  dependent  on  wind  direction.  If  the  wind  directions  are  not  right  on  the 
line  from  the  sources  to  the  receptors,  then  the  receptors  may  not  be  "seeing" 
the  plumes  from  the  sources.  The  diffusion  coefficients  (sigma  y and  sigma  z) 
are  supposed  to  handle  some  of  this  fanning,  but  may  not  be  large  enough  to 
account  for  what  is  seen  in  Billings. 

The  Department  did  try  to  validate  a model  so  that  the  results  compared 
better  with  the  actual  concentrations.  This  involved  changing  some  of  the  input 
parameter  for  the  sources.  For  instance,  since  effective  plume  height  seemed 
too  high,  the  Department  lowered  stack  exit  temperatures  to  293°K  so  that  no 
plume  rise  occurred.  In  some  cases,  the  physical  stack  heights  were  also 
lowered.  The  Department  also  raised  emission  rates  at  Montana  Power  and  Exxon 
to  reflect  a balance  of  concentrations  seen  from  those  directions.  The 
regressions  (comparing  actual  to  modeled  concentrations)  improved  slightly  but 
the  juggling  of  figures  seemed  inappropriate. 

The  Department  concludes  that  Complex  I and  II  are  not  able  to  model 
sources  in  Billings  with  its  complex  terrain  and  meteorological  features.  This 
conclusion  also  applies  to  any  Gaussian  dispersion  model. 
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XX.  SOURCE  CONTRIBUTION 


This  section  is  intended  to  show  how  the  Department  calculated  source 
contribution  by  the  Billings  and  Laurel  industries  on  the  annual  average  SO^ 
concentrations  at  Coburn  Road  and  Laurel -BN  sites  during  the  study  year. 

In  order  to  determine  source  contribution  from  the  six  industries  on  the 
annual  average  at  Coburn  Road,  a couple  of  procedures  were  applied.  First,  in 
order  to  get  an  idea  of  the  magnitude  of  each  source,  the  Department  calculated 
percentages  based  on  mass  emissions.  (It  is  noted  that  this  procedure  considers 
no  knowledge  of  meteorology,  dispersion  and  location  of  receptors  with  respect 
to  the  location  of  the  sources  and  assumes  that  the  concentrations  are  directly 
proportional  to  emissions.)  The  percentages  listed  on  Table  148  for  the  five 
Billings  industries  were  used  for  this  determination.  In  order  to  figure  in  the 
contribution  by  Cenex,  the  Department  arbitrarily  assigned  that  plant  a five 
percent  contribution.  The  Billings  plants  were  then  multiplied  by  95  percent 
(or  0.95).  The  results  of  this  determination  are  shown  below: 
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A second  procedure  was  also  used  to  calculate  source  contribution  on  Coburn 
Road.  This  procedure  used  the  direction  analysis  results  to  proportion  out 
subareas  (A,  B and  C)  and  then  used  mass  emissions  to  proportion  out  individual 
sources  within  those  subareas. 
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This  second  procedure  followed  some  individual  steps.  The  first  step  is  to 
add  up  the  percent  strength  factors  for  the  SW  directions  and  for  the  NE  direc- 
tions. Table  156  shows  the  summaries  of  the  four  study  months.  The  percent 
strength  factors  from  the  SW  direction  show  the  contribution  from  subareas  A and 
B while  the  percent  strength  factors  from  the  NE  directions  show  the 
contribution  from  subarea  C. 

The  next  step  is  to  average  the  SW  percent  strength  factors  and  to  average 
the  NE  percent  strength  factors  for  the  four  study  months.  The  Department 
assumes  that  the  averages  of  the  four  study  months  are  representati ve  of  the 
annual  average.  A summary  of  this  averaging  is  as  follows: 


SW 

NE 

SW  + NE 

January 

56.7% 

38.7% 

95.4% 

May 

51.0 

33.4 

84.4 

August 

34.6 

48.9 

83.5 

November 

81.4 

9.7 

91.1 

Average 

55.9 

32.6 

88.5 

The  average  percent  strength  factors  for  the  SW  directions  equaled  55.9% 
and  the  average  for  the  NE  direction  equaled  32.6%.  The  total  of  the  two 
equaled  88.5%  or  explained  approximately  90%  of  the  SO,-,  seen  at  the  Coburn  Road 
site. 

The  third  step  is  to  set  up  a ratio  and  equate  it  with  the  sources  in  those 
areas.  The  ratio  and  the  equation  is  as  follows: 

SW  = 55.9  = 1.71  = A + B (1) 

NE  ITT  ~~~Z 

The  fourth  step  is  to  equate  the  sum  of  SO^  from  each  of  the  subareas  with 
total  annual  average  at  Coburn  Road.  This  is  as  follows: 

A + B + C = .029  ppm  (2) 
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A situation  presents  itself  by  having  two  formulas  and  three  unknowns:  If 

one  of  the  unknowns  can  be  determined  independently  of  the  others  and  substi- 
tuted into  equations  (1)  and  (2)  above,  then  the  other  two  values  can  be  solved 
using  simultaneous  equations.  This  leads  to  the  next  step. 

The  fifth  step  is  to  solve  for  A:  The  amount  or  fraction  of  S0o  from 

L. 

subarea  A (The  Laurel  Cenex  plant)  that  contributes  to  the  annual  average  at 
Coburn  Road  for  this  study  year.  This  is  done  by  making  two  assumptions.  These 
are  as  follows: 

1.  All  of  the  SO^  from  the  SW  sector  at  Taft  School  is  from  subarea  A. 

2.  The  amount  of  SO^  at  Taft  School  from  subarea  A is  equal  to  the  amount 
of  SO^  at  Coburn  Road  from  subarea  A. 

The  fraction  of  SQ^  from  subarea  A at  Taft  School  is  found  by  adding  up  the 
percent  strength  factors  from  the  SW  directions.  The  concentration  value  for 
the  SW  sector  is  found  by  multiplying  the  sum  of  the  percent  strength  factors  by 
the  monthly  average.  This  value  is  divided  by  the  monthly  average  at  Coburn  to 
show  percent  contribution  from  subarea  A at  that  site.  The  figures  for  the  four 
study  months  are  averaged  to  give  approximate  annual  averages.  Tables  157  to 


160  show 

these  results  and  a 

brief  summary 

is  presented  below: 

Taft  SW 

Taft 

Taft 

Percent 

Percent 

Monthly 

SW 

Taft  (SW)/ 

Strength  Factor 

Average 

Concentration 

Coburn  (total) 

January 

20.5% 

21.8  ppb 

4.4  ppb 

8.6% 

May 

16.2 

6.1 

1.0 

15.6 

August 

29.2 

12.4 

3.6 

13.0 

November 

56.5 

12.2 

6.9 

12.8 

Average 

30.6 

13.1 

4.0 

12.5 

From  this  step  it  is  concluded  that  the  value  for  A (the  amount  of  S02  from 
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Laurel  at  Coburn  Road)  is  equal  to  .004  ppm  (or  4 ppb).  This  step  also  shows 
that  12.5%  of  the  SO^  at  Coburn  Road  is  from  the  Cenex  plant. 

With  A = .004,  equations  (1)  and  (2)  reduce  to  the  following: 

1.71  = .004  + B 

T 

and 

.004  + B + C = .029 

Solving  the  two  equations  simultaneously  gives  B = .0143  or  approximately 
.014  and  C = .0107  or  approximately  .011.  In  summary,  the  contributions  from 
the  three  subareas  are  as  follows: 

A = .004 
B = .014 
C = .011 

The  sixth  and  final  step  is  to  use  the  partial  contributions  by  subarea  and 
the  mass  emissions  presented  on  Table  147  to  calculate  source  contribution.  The 
results  are  as  follows: 

MPC  = .014  X (5460/9626)  / 100  = .00794  = 27.4% 

CON  = .014  X (3603/9626)  / 100  = .00523  = 18.0 

GTW  = .014  X ( 563/9626)  / 100  = .00081  = 2.8 

EXX  = .011  X (8269/10796)  / 100  = .00842  = 29.0 

MTS  = .011  X (2527/10796)  / 100  = .00257  = 8.8 

CEN  = .004  X (8794/8794)  / 100  = .00400  = 13.8 

.02897  99.8 

The  following  table  presents  the  results  of  the  two  estimation  methods. 
Since  there  is  a degree  of  uncertainty  in  the  determinations,  the  Department 
shows  a range  of  values  in  which  the  actual  source  contributions  would  probably 
occur.  The  table  also  presents  single  percent  numbers  that  can  be  understood  as 
the  Department's  best  estimate  of  the  concentrations  by  the  six  industries  on 
the  annual  average  at  Coburn  Road  during  the  study  year. 
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Estimation 
Method  1 

Estimation 
Method  2 

Range 

Best 

Estimate 

Montana  Power 

25.3% 

27.4% 

25“ 33% 

31% 

Conoco 

16.7 

18.0 

13-19 

14 

Great  Western 

2.5 

2.8 

2-4 

3 

Exxon 

38.2 

29.0 

26-34 

30 

Montana  Sulphur 

11.6 

8.8 

8-12 

10 

Cenex 

5.0 

13.8 

10-14 

12 

99.3 

99.8 

84-116 

100 

The  Department  also  wished  to  calculate  the  source  contributions  on  the 
Laurel-BN  site  in  Laurel.  This  was  done  by  a method  similar  to  what  was  used 
for  the  Coburn  Road  site. 

The  first  step  is  to  add  up  the  percent  strength  factors  for  the  SW 
directions  and  for  the  NE  directions  for  each  of  the  four  study  months.  The 
percent  strength  factors  from  the  SW  directions  show  the  contribution  from 
subarea  A while  the  percent  strength  factors  from  the  NE  directions  show  the 
contribution  from  subareas  B and  C.  Table  161  shows  these  results. 

The  second  step  is  to  average  the  SW  percent  strength  factors  and  to 
average  the  NE  percent  strength  factors.  The  Department  assumed  that  the 
average  of  the  four  study  months  are  representative  of  the  whole  study  year.  A 
summary  is  presented  as  follows: 


37 


SW 

NE 

SW  + NE 

January 

58.3% 

38.1% 

96.4% 

May 

65.5 

21.3 

86.8 

August 

59.8 

19.2 

79.0 

November 

81.7 

9.6 

91.3 

66.3 

22.0 

88.3 

The  average  percent 

strength  factors  for  the  SW 

directions  equaled  66.3% 

and  the  average  for  the 

NE  directions 

equaled  22.0%. 

The  total  of  the  two  gave 

88.3%  or  explained  approximately  90%  of  the  SOo  seen  at  the  Laurel-BN  site. 

For  the  third  step,  the  Department  assumed  that  all  of  the  SO^  from  the  NE 
direction  came  from  the  Billings  industries  (or  22%  of  the  annual  average  at 
Laurel-BN  site).  By  using  the  22%  figure,  the  sources  in  Billings  can  be 
proportioned  by  mass  emissions  or  using  percentages  presented  in  Table  148.  The 
following  summarizes  these  results. 
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Since  there  is  a degree  of  uncertainty  in  the  determinations,  the  Depart- 
ment figured  a range  of  values  in  which  the  actual  source  contribution  would 
probably  occur.  The  following  table  shows  this  range  and  the  Department's  best 
estimate  of  the  contributions  by  the  six  industries  on  the  annual  average  at  the 
Laurel-BN  site  during  the  study  year. 
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Estimation 

Method 

Range 

Best 

Estimate 

Montana  Power 

5.8% 

5-9% 

6% 

Conoco 

3.8 

3-7 

4 

Great  Western 

0.6 

0-3 

1 

Exxon 

8.9 

8-12 

9 

Montana  Sulphur 

2.7 

2-6 

3 

Cenex 

78.0 

70-80 

77 

99.8 

88-117 

100 
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XXI.  DISCUSSION 


The  purpose  of  this  section  is  to  briefly  discuss  some  of  the  aspects  of 
the  report  and  to  emphasize  some  points  that  the  Department  feels  are  important 
to  the  final  determination. 

The  objectives  of  the  stipulation  were  to  (1)  determine  air  quality  S0? 
levels,  (2)  validate  a model  and  (3)  use  the  data  and  model  to  determine  source 

contributions. 

The  air  quality  levels  during  this  study  period  showed  that  there  were  no 
violations  of  any  Federal  SO^  standards.  However,  there  were  violations  of  the 
Montana  24-hour  and  annual  SCL  standards.  The  Coburn  Road  site  had  the  highest 
annual  average  in  Billings  (.029  ppm)  and  the  Laurel-BN  site  had  the  highest 
annual  average  in  Laurel  (.010  ppm).  The  Montana  annual  standard  is  set  at  .020 
ppm. 

The  second  objective  (to  validate  a model)  was  perhaps  simply  and  optimis- 
tically stated.  It  was  not  exactly  clear  in  the  stipulation  what  was  meant  by 
the  word  "model."  However,  as  a first  approach,  the  intent  did  seem  to  mean 
Gaussian  plume  dispersion  models.  In  fact,  emission  and  meteorological  data  was 
collected  which  would  be  input  to  such  models.  The  Department  endeavored  to 
validate  Gaussian  plume  dispersion  models  but  the  effort  was  unsuccessful  and 
the  models  were  invalidated  (or  disproved)  for  the  Billings  and  Laurel  areas. 

Since  Gaussian  plume  dispersion  models  were  disproved  for  the  area,  the 
stipulation  did  seem  to  allow  the  Department  to  apply  alternative  methods 
(models)  to  determine  source  contribution.  The  Department  used  direction 
analysis  and  mass  emissions  to  calculate  source  contribution.  This  was  shown  in 
the  previous  section. 
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Although  the  models  were  disproved  for  the  Billings  area,  it  should  not  be 
assumed  that  all  modeling  endeavors  give  poor  results.  Gaussian  plume  dispers- 
ion models  are  still  the  best  tool  available  to  show  possible  impacts  by 
proposed  sources.  The  Department  feels  that  they  should  definitely  be  retained 
as  part  of  the  permitting  process.  The  authors  of  the  models  maintain  that  they 
do  not  include  calculations  for  the  transitional  phenomenon  of  fumigation  or  the 
breakup  of  inversions.  It  is  no  wonder  that  the  models  were  unable  to  predict 
concentrations  from  the  Billings  and  Laurel  industries. 

The  Department  notes:  If  the  same  sources  existed  on  a flat  terrain  with 
the  same  set  of  emissions,  then  there  probably  wouldn't  be  any  violations  of  the 
ambient  standards.  This  is  essentially  what  the  models  are  predicting.  Given 
the  intervening  terrain  and  meteorology,  a situation  creates  itself  whereby  SO^ 
concentrations  build  up  or  impact  the  ground  and  cause  violation  of  standards. 

In  order  to  give  more  understanding  of  the  inversion  breakup  phenomenon. 
Figure  76  shows  a sequence  of  valley  cross  sections  illustrating  the  breakup  of 
a valley  temperature  inversion  beginning  at  sunrise  and  ending  at  the  time  of 
inversion  destruction.  This  figure  is  depictive  of  what  happens  in  the  Billings 
and  Laurel  areas  during  all  parts  of  the  year.  (This  figure  is  referenced  to  a 
paper  written  by  C.  David  Whiteman  and  K.  Jerry  All  wine  which  is  entitled  " A 
New  Pollutant  Dispersion  Model  for  Application  in  Mountain  Valleys.")  The 
inversion  is  destroyed  by  the  upward  growth  from  the  ground  of  a convective 
boundary  layer  and  the  subsidence  or  sinking  of  the  temperature  inversion  layer 
after  sunrise. 

The  stipulation  addressed  the  possibility  of  raising  stack  heights  to 
control  emissions  or  to  lower  air  quality  levels.  Modeling  assumptions  also 
base  their  predictions  on  stack  heights  and  plume  rise.  It  is  the  Department's 
opinion  that  stack  heights  in  Billings  do  little  to  control  or  lower  ambient 
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concentrations.  For  example,  plume  impaction  by  subarea  C (Exxon  and  Montana 
Sulphur)  is  seen  at  the  Lockwood  Park  site  which  is  only  a half  a mile  away. 
Also,  the  Montana  Power  Corette  stack  impacts  the  Coburn  Road  site  which  is 
about  a mile  away.  Again,  stack  heights  seem  to  do  little  to  lower  ambient 
concentrations  in  Billings  and  raising  stack  heights  should  not  be  an  acceptable 
method  for  controlling  S0?  levels. 

The  signees  of  the  stipulation  were  also  interested  in  whether  there  were 
any  areas  that  would  show  higher  levels  than  those  measured  by  the  ambient 
monitors.  This  would  have  been  accomplished  by  adding  additional  receptor 
locations  in  the  models.  Since  the  models  did, not  work,  the  Department  makes 
the  decision  based  on  weather  patterns  and  observed  ambient  levels. 

The  Department  feels  that  the  Coburn  Road  site  is  probably  the  worst  case 
site  in  the  Billings  and  Laurel  area.  The  Coburn  Road  site  receives  SO^ 
emissions  from  both  the  NE  and  5W  directions  and  also  from  above  during  inver- 
sion breakups.  The  Coburn  Road  site  is  also  an  elevated  site  (sits  high  on  a 
ridge)  and  is  at  a height  where  buildups  occur  during  calm  weather  conditions 
(close  to  the  inversion  level).; 

The  Department  wishes  to  point  out  that  there  are  areas  other  than  the 
immediate  Bi 1 1 ings-Laurel  area  that  are  affected  by  the  emissions  from  the 
Billings  and  Laurel  industries.  It  is  the  Department's  judgement  that  the  SC^ 
may  stretch  along  a line  over  the  Yellowstone  River  to  points  as  far  away  as 
Columbus  to  the  SW  and  Pompey's  Pillar  to  the  NE.  The  concentrations  judged  at 
these  distant  locations  could  be  as  high  as  .005  ppm  on  an  annual  average. 

The  Department  also  feels  that  there  are  no  significant  SO^  sources  in  the 
Billings  area  other  than  the  industries  that  this  report  addresses:  Montana 

Power,  Conoco,  Great  Western,  Exxon,  Montana  Sulphur  and  Cenex.  There  may  be 
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some  small  miscellaneous  sources , but  these  would  probably  contribute  less  than 
a half  of  percent  of  all  the  emissions. 
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XXII.  SUMMARY  AND  CONCLUSIONS 

1.  The  ambient  air  quality  data  collected  during  the  study  year  (October  1, 


1981 

to  September  30,  1982)  shows 

violations  of  state  annual  and  twenty-four 

hour 

standards.  There  were  no  violations  of  any  Federal  standards. 

2. 

The  source  contribution  for 

the  Coburn  Road  site  (annual  average  equal  to 

.029 

ppm  SO^)  is  as  fol lows: 

Montana  Power 

31% 

Conoco 

14 

Great  Western 

3 

Exxon 

30 

Montana  Sulphur 

10 

Cenex 

12 

100% 

3. 

The  source  contribution  for 

the  Laurel-BN  site  (annual  average  equal  to 

.010 

ppm  SO^)  is  as  follows: 

Montana  Power 

6% 

Conoco 

4 

Great  Western 

1 

Exxon 

9 

Montana  Sulphur 

3 

Cenex 

77 

100% 

4.  Gaussian  plume  models  were  proved  invalid  for  the  Billings  and  Laurel  area 
and  cannot  be  used  to  determine  source  contribution. 
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5.  Federal  annual  SOg  violations  might  occur  if  the  Montana  Power  Corette 
plant  ran  for  the  whole  year  or  if  there  were  some  increases  in  production  at 
the  oil  refineries. 
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100  copies  of  this  public  document  were  published  at  an 
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which  includes  $325.00  for  printing  and  $.00  for  distribution. 


